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AN ASYMPTOTIC INVESTIGATION 

OF THE STATIONARY MODES OF INSTABILITY OF THE BOUNDARY LAYER 

ON A ROTATING DISC 
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Abstract 

We investigate high Reynolds number stationary instabilities in the 
boundary layer on a rotating disc. The investigation demonstrates that in 
addition to the inviscid mode found by Gregory, Stuart, and Walker (1955) 
at high Reynolds numbers, there is a stationary short wavelength mode. 
This mode has its structure fixed by a balance between viscous and Coriolis 
forces and cannot be described by an inviscid theory. The asymptotic 
structure of the wavenumber and orientation of this mode is obtained, and a 
similar analysis is given for the inviscid mode. The expansion procedure 
provides the capacity of taking non-parallel effects into account in a 
self-consistent manner. The results are compared to numerical calculations 
and experimental observations. 
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1. Introduction 


In recent years there has been much interest m the manner m which 
three-dimensional boundary layers become unstable. Much of this work has 
been motivated by the need to understand the instability mechanisms which 
are operational in the boundary layer on a swept wing. This research has 
been directed towards the development of laminar flow airfoils and the 
possible instability of the flow to Gortler vortices and crossflow 
vortices, while taking Tollmien-Schlichting waves into account. Thus, Hall 
(1985) considered the Gortler vortex instability in a weakly three- 
dimensional boundary layer and found that an asymptotically small spanwise 
velocity component is sufficient to prevent the Gortler mechanism occurring 
at finite Gortler numbers. 

The crossflow mechanism has also been the subject of many 

investigations (see, for example, Gregory, Stuart, and Walker (1955), 
Cebeci and Stewartson (1980), Malik, Wilkinson, and Orszag (1981), and Reed 
(1985)). This instability mechanism occurs when the effective velocity 
profile in an Orr-Sommerf eld approximation to the linear instability 
equations has an inflection point where the velocity field vanishes. The 
importance of such a profile was explained by the inviscid analysis of 
Gregory, Stuart, and Walker (hereafter referred to as GSW) in the context 
of the rotating disc problem. 

The noteworthy feature of this profile is that it can support a 

stationary vortex pattern relative to the disc. GSW showed that the normal 
to the vortex boundaries made an angle <(> of about 13° to the radius 
vector. This was found to be in excellent agreement with their 
experimental observations, but the number of vortices predicted by the 
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theory was found to be too large by a factor of about 4. The latter 
discrepancy has been attributed to viscous effects, but the reason why the 
angle <j> should not also be significantly altered by such effects is not 
clear. The asymptotic investigation of the GSW mode which we will give in 
Section 3 will shed light on this question. 

A recent parallel flow numerical investigation by Malik (1983) found 
that the point at infinity of GSW in the wavenumber-Reynolds number plane 
is connected to a curve corresponding to stationary modes at finite 
Reynolds number. However, the angle c(> varies along the curve and the 
critical Reynolds number corresponds to <|> ~ 11°, and there is also a lower 
branch on which asymptotes to about 39° when the Reynolds number is 
larger. 

The first purpose of this paper is to set up a rational framework which 
can take non-parallel effects into account at large Reynolds numbers. The 
second aim is to provide an analytical method of producing the wavenumber- 
Reynolds numbers dependence of the upper and lower branch modes. Since our 
analysis is applicable to any three-dimensional boundary layer, our 
calculations enable the likely stationary vortex patterns in such flows at 
high Reynolds numbers to be predicted analytically. We will see that the 
lower branch mode corresponds to the case when the effective velocity 
profile has zero shear stress at the wall and the disturbance takes on a 
triple deck structure. The development of an asymptotic theory will also 
enable nonlinear effects to be investigated in a self-consistent manner. 
Such an investigation is beyond the scope of the present paper but is 
clearly necessary in order to explain why the upper branch mode is 
apparently almost always the only one to be observed experimentally. The 
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asymptotic theory of the lower branch mode is also relevant to short 
wavelength instabilities of Stokes layers. The procedure adopted in the 
rest of the paper is as follows: m Section 2 we formulate the instability 
equations; in Sections 3 and 4 we develop asymptotic theories for the upper 
and lower branch modes. Finally, in Section 3 we draw some conclusions. 


2. Formulation of the Problem 

We consider the flow of a viscous fluid of kinematic viscosity v in 
the region z > 0. The motion of the fluid is induced by the steady 

rotation with angular velocity 0, of the plane z = 0 about the z axis. 
We take cylindrical polar coordinate (r, 0, z) with r and z having 
been made dimensionless with respect to some reference length l. The 
Reynolds number R for the flow is defined by 



( 2 . 1 ) 


and if the axes rotate with the plane, then the basic steady velocity field 
is 

£ = it, = £ft(ru(Rz ^ ) , rv(Rz ^ ) , R ^ 2 w(Rz ^ ) ). (2.2) 


Here the functions u, v, and w are determined by 

u 2 - (v + l) 2 + u" w - u" = 0, 
2u(v + 1) + v' w - v" = 0, 


(2.3a) 


(2.3b) 
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2u + w' = 0, 


(2.3c) 


where the prime denotes differentiation with respect to z. The 
appropriate boundary conditions are 


u = 0, v = 0, w = 0, z = 0 
u -► 0 , v -*• -1 , z 


(2.4) 


We now perturb the above flow by writing 


u = Ug + ft£U((r,0,z), V(r } 9,z), W(r,0,z)) 


(2.5) 


where U, V, and W are small and steady. The expression (2.5) is then 
substituted into the Navier-Stokes equations in the rotating frame and 
linearized to give 


r — 3 

|ru 


h * v I? + R ~ '' 2 " lit" + uU - + rW if 


V? — 3 


3u 




2_ 9V 

2 ae 




(2.6a) 


— 8 - - 9 -*- R 2/ 2 w |_Jv + uV + 2{v+l }U + rW |^- 


{-f7 + -flT + 


9z 


I II 

r 90 


+ 


I J/.V + — — - — } 

R 2 90 2 > 


> 


(2.6b) 
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{ru |- + v + R _ w -| — }W + R~ W 4^ 
1 3r 30 3z J 3z 




(2.6c) 


where 



1 

2 

r 



+ 


±f- + 

r 3r 



> 


and P is the nondime ns tonal pressure perturbation. The equation of 
continuity then becomes 


3£ + U + J_3V_ + 3W 
3r r r 39 3z 


0 . 


(2.7) 


Finally, we must solve (2.6), (2.7) subject to the no-slip condition at the 
wall, whilst sufficiently far away from the wall we insist that the 
disturbance decays to zero. However, we shall see below that the length 
scale for this decay to zero will depend on the type of disturbance under 
consideration. 


3. The Inviscid Modes 

From the inviscid theory of GSW we expect these modes to have 
wavelengths scaled on the boundary layer thickness. Thus, we must consider 
modes with a length scale of order R ^2 xn the r and 0 directions. 
It is convenient to define the small parameter e by 


e = 


r - 1/6 
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and then we write 

r 

U = u(z) exp iy { j a(r,e)dr + 93(e)} (3.1) 

e 

together with similar expressions for V, W, and the pressure perturbation 
P. The wavenumber a will, m general, be complex and is determined in 
terms of e and 3. However, we will restrict our attention to neutral 
disturbances and find a and 3 such that the flow is neutrally stable at 
the position r. We now expand a, 3 as 

a = (Xq + ecij + >•••> (3.2a) 

3 = 3q + eBj + »••• . (3.2b) 

The disturbance structure in the z direction is then fixed by the following 

considerations. Firstly, from the results of GSW we anticipate that there 

3 

will be an inviscid zone of depth 0(e ). In order to satisfy the no-slip 

condition on the velocity at the wall, a viscous layer must exist. The 

4 

thickness of this layer is then found to be 0(e ) by balancing the 
convection and diffusion terms in the disturbance equations. 

In the inviscid zone we expand u, v, w, and p m the form 


u = u q ( 4) + euj(c) + ,*•*, (3.3a) 
v = v Q (s) + ev^s) + ,•••, (3.3b) 
w = w Q (?) + ew^s) + ,•••, (3.3c) 
P = Pq(s) + E Pi(?) + »•**» (3.3d) 



-7- 


where 5 = ze . The above expansions are then substituted into (2.6), 
(2.7) with replaced by + ij {a Q + eo^ + ,*•*,} and with -jL- 

i r 1 E -4 

replaced by — j {g^ + eg^ + ,•••,]■- If we equate terms of order e , we 

e 

obtain 


1 UUQ + rw 0 u' = - 1« 0 p Q 

= — , lB 0 

luv o + ™o v = ' ~r p 0 


(3.4a) 

(3.4b) 


luw 0 = “ p 0 


i6v c 

ict 0 u 0 + ~ 


+ w' = 0 , 


(3.4c) 

(3.4d) 


where u = cIq ur + g^ v. If we eliminate uq, Vq and the pressure from 
the above equations, we find that wq satisfies 


U K " Y o w q] ■ u " w o = 0 


(3.5) 


with u now acting as the 'effective' or 'equivalent' two-dimensional 

s 2 

2 2 p 0 

velocity profile whilst ^0 = a 0 + ~ t * ie e ff ectlve wavenumber. Thus, 
wq satisifies Rayleigh's equation and is determined as an eigenvalue 

when (3.5) is solved subject to 


Wq = 0, ? = 0,». 


(3.6) 


We further note that ®q/3q is chosen such that u and u" vanish at the 
same nonzero value of £ = £; in this case (3.5) has no singularity at 
£ = £. The eigenvalue problem was solved using central differences; we 
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obtained 


= 1.16, 


(3.7a) 


(3.7b) 


= 1.46. 


The eigenfunction wq normalized with w' = 1 at £ = 0 is shown in 
Figure 1. 

Having calculated Wq we can use (2.8), (2.9) to solve for uq, vq, 
and pq; however, it suffices here to say that when 5 -*■ 0 


6 0 V Q i 


1 {“o u o + “ — } + “ w o (0) - 


(3.8) 


Before proceeding to the next order in the inviscid zone, we calculate the 
zeroth-order solution in the wall layer. If we write 


K = e 4 z, 


then in the wall layer u expands as 


u = e u Q % + ,•••, 


and v, w are expanded m a similar manner. The disturbance velocity and 


pressure now are written as 
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u = U 0 (5) + eU i ^ + > ' * * » 
v = V Q (C) + eV^O + 

W = £Wq(§) + £ 2 W^C?) + ,•••, 

s 

p = £Pq( 5) + e 2 P^?) + • 


( 3 . 9 ) 



Wq(0) A'(0) 

W ~ w' £ + 

0 0^ 00 » 

Y / A (s)ds 

0 1 

so that the order e inviscid zone normal velocity component must 
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satisf y 


w'(0)A'(0) 

w ->• + , £ * 0. (3.13) 

1 00 

Y / A (s)ds 

0 1 

We now turn to the next order problem in the inviscid zone. Thus, we now 
-3 

equate terms 0(e ) in the inviscid zone disturbance equations, and we 

obtain a set of equations similar to (3.4) with (uq, Vq, Wq, Pq) replaced 
by (uj, Vj , w^ , p^) but now having inhomogeneous terms. If we repeat the 
manipulations carried out on (3.4) m order to get (3.5), we obtain 


S 0 0 1 ( 


u K _ Y o W J " u " w i = 2u a o a i + — rr o 


+ a, - 


8 1 “0 


1 0. 


r< u" - ^ 


w 0 . (3.14) 


The second term on the right-hand side of (3.14) causes w^ to have a 
logarithmic singularly at £ = £; this can be removed in the usual way by 
introducing a critical layer at £ = £. We can formally write down a 
solution of (3.14) which satisfies w^(“) =0 in the form 


0n 0 


W 1 = 2 “ 0 °1 + 


V 1 w o (?) / / w O (0)d0 

r 1 w Q (n)J 


+ a, - 


a o e i| 


1 0C 


rw Q (4) 


f \~)f 

= w 0 (ti) » 


( 0 ) 


u"(0)u(0) - u(0)u(6)d0 


=2 

U 


(3.15) 


when £ is a constant with £ > £. The above solution is valid for £ < £ 
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lf the path of integration is deformed appropriately into the complex plane 
near £ = It can then be shown from (3.15) that 


^ J ^ e l B 0 ) X 1 ^(“l_ a O e i( r T ,, 

W 1 ( 0) 2 j a o a i + ~ _ ^r~j^T(5T I 2’ (3, } 


where 


h ■/ ”0 (9)d6 


(3.17a) 


„ f 2 u" u - u" u 

= 8 0 J W 0 =2 


(3.17b) 


where the path of integration m (3.17b) is deformed below the singularity 
at 5 = 5". The matching condition (3.13) produces the eigenrelation, 


A iq( w 0 (0)2 ) 

00 

Y J A.(s)ds 
0 1 


= 2 a Q a L + 




( 6 0 


a o B i) T 

2 ( r l 2‘ 


(3.18) 


Our calculations showed that 


I 1 = .094, I = .058 - .029i, 


and using the well known values for A'q, f A^(s)ds, we obtained 


a 0 + 


e 0 e i ,, -1/3 

T~ = ~ 14 * r V 


!i - r - 2V • r- 1/3 . 


B ° »0 f 
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The above equations can be solved for , g^; however, it is more useful 
to evaluate 


2 2 2 

/a + r/r = y q + 


a 0 a l 


g Q 



= 1.16 


14.4 -1/3 A 

?76 r + "" 


(3.19) 


which we interpret as the 'effective' wavenumber of the disturbance. 
We now define the wave angle <f) by 



^o + K 

e o r ( 4 * * * 8 o 


6 1 a o, 


e + 


4.26 + 


29. r 


-1/3 


, 1/6 


(3.20) 


Thus, we have calculated the first correction terras to the classical 
results of GSW. The sign of the correction term in (3.20) has some 
important consequences which we will discuss in Section 5. 


4. The Wall Modes 

We have seen in the previous section that the 'effective' velocity 

profile for a three-dimensional disturbance with wavenumbers a and g m 
the r and 0 directions is aur + gv. The inviscid modes are such that 

our + gv and au"r + gv" vanish simultaneously. It is easy to show that 

lower branch disturbances having a triple deck structure of the type 

discussed by Smith (1978) for Blasius flow can also exist. However, such 
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modes are necessarily time-dependent with a, 8 real if the effective wall 
shear au'r + f3v' does not vanish. Therefore, we choose to look for 
stationary modes for which the effective wall shear vanishes at zeroth 
order. 

It is easy to show that the appropriate triple-deck structure is based 
on the small parameter e now defined by 


e 




( 4 . 1 ) 


and the lower, main, and upper decks are of thickness e , e , and e 
respectively. The disturbances structure m the main and upper decks is 
essentially the same as that found by Smith (1978) who investigated lower 

branch disturbances to Blasius flow. The wavenumbers in the r and 0 

-4 

directions are now 0(e ); we therefore write 

r 

U = U(z)exp ^ { / a(r,e)dr + 98(e)}, 
e 

together with similar expressions for V, W, and P. We define £, £, Z by 


C 


= 

9 * 
e 




The wavenumbers then expand as 


a 


+ 


2 

e 



9 


• • • 


9 


8 = 8 0 + e 


+ e 3 B 2 + 


(4.2a,b,c) 


(4.3a) 


(4.3b) 
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Here we have anticipated that the order e terms are zero, and we again 
seek a^, g^, etc. such that the flow is neutrally stable at the location 
r. In the upper deck u = 0, v = -1 , and U expands as 


U = e ' 


V 2) 


+ e 


Uj(z) 


and V, W, and P have similar expansions. We found that the zeroth- 
order equations to be solved in the upper deck are 


^0 P 0 _ dP 0 

5 0 U 0 " a 0 P 0’ 3 0 V 0 r ’ l6 0 W 0 " dz ’ 


6 0 v o dW o n 

la 0 U 0 + 1 r + d T = °’ 


and the solution of this system which decays to zero when z -»■ is 


P 0 = Ce 


-Y 0 z 


U o = ^ Ce 


a o „ -Yo 2 


(4.4a,b) 


V n = — e 
0 r 


C “*0 2 


1Y o „ " Y o 2 


w o=V Ce 


(4.4c,d) 


where 


/ 2 7 2 . 2 

Y 0 /a Q + g Q /r. , 


and C is an unknown function of r. 


In the main deck the disturbance expands as 
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U = j u Q (?) + u^c) + ,•••, 

V = J v q(c) + v j(?) + »•••» 

3 2 

W=e w Q (?)+e Wj (5) +,•••, 
P = e 3 C + e 2 P^c) + 


where we have anticipated that P is independent of £ to order e and 

therefore equal to C. Substituting into (2.6), (2.7) we find that uq, 

Vq, wq satisfy 


ioiq ru u Q + i 8 q v Uq + ru' Wq = 0 , 
iccq ru Vq + i 0 q v Vq + rv' w Q = 0, 


“0 u 0 + r 


1 B 0 

— v 0 + dT 


0, 


and the solution of this 
is 

u 0 


v 


0 


system which matches with the upper deck solution 


Cr ^0 u ' 


Cr ^0 v ' 


(4.5a) 

(4.5b) 


w 


0 


—■ (a 0 ru + Sq v). 

e o 


(4.5c) 


We note from (4.5c) that 
when £ -*■ 0; however, unless 


Wq in fact satisfies the no-slip condition 
u' and v" both vanish at 5 = 0 the other 
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velocity components are nonzero there. If we choose cxq and such 


that 


— — 
a Q u'( 0 ) + -f- v'( 0 ) = 0 , 


(4.6) 


a 0 r 


*0 v 0 


which gives — = 1.207, then c*q + — — — + 0 when 5 ->■ 0. It is 

the imposition of the constraint (4.6) on the effective velocity profile 
which enables us to find stationary disturbances. If we expand u, v for 
small £ and write 5 = £/e, we have 


— — 2—2 — 3 

u = £Uq £ + e UjC +u 2 ? 


(4.7a) 


v = ev Q C + e 2 Vj ? 2 + v 2 £ 3 + ,••*, 


(4.7b) 


— U J ( 0 ) * — v J ( 0 ) 

when u , = — ) , v. . = i — ; for i = 1,2,»»* . In order to match 

J-l j! J-l j! 

with the solution (4.5), written in terms of £ using (4.7), we therefore 
expand the lower deck disturbance in the form 


r Yn c 


U = 


eB^ 


u Q + 2 eu x £ + 


U ,(€) 

+ — + u q(?) + eU^(?) + ,•••, (4.8a) 


V = 


^0 


ee; 


v Q + 2ev l 5 + 


V_j(5) 


+ + eV^(£) + (4.8b) 


iYn e C 


W = - 


(“0 U 1 r + 8 o V 1 )?2 + 


+ e Wj (?)+,•••, 


(4.8c) 


P = e P^S) + ••• • 


(4.8d) 
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We must now substitute the above expansions into the disturbance 
equations and solve for (U_j , V..^), (Uq, Vg), (U^ , , Wj , Pj), etc. From 
the continuity equation we obtain immediately that 


-1 


I2-U 

6 0 r -1’ 


= - 


U n » 


0 Bg r “0 


(4.9a) 

(4.9b) 


0 

where - — satisfies (4.6). From the radial momentum equation we obtain 
6 0 


- 2 d2 °-l 

- i[rdg u 1 + Bg vjs U_J + — = 0, (4.10a) 

d? 

d 2 U 

i[ra 0 u 1 + Bg vj? U Q + 2 = " ru 0 W 1 + t ra 0 U 2 + 6 0 V 2^ U 0’ 

d? 

(4.10b) 


which must be solved subject to 


-1 


^0 Cu o 


, U Q = 0, % = 0, 


(4.11) 


u _!> % > 0, 5 + 


The function U_ 1 is given by 


-1 


u o Y o Cr ufo, n a 1/4 c) 

2 U(0,0) 

Pn 


(4.12) 


where 
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A = i{o Q ru x + B 0 Vj} 


(4.13) 


and U(0, /2 A is a parabolic cylinder function. The functions Uq, 

Vq cannot be determined unil is calculated, the latter function can 

be found by considering the next order approximation to the radial and 


azimuthal momentum equations. If we multiply these equations by ioiq and 


iB 


0 


respectively and add them we obtain: 


i<a 


“ 2 “l . S 0 “ 2 “l / . 2 




r ~~2 J + Y 0 P Q + 2la 0 V -1 r 


2 i8 0 ^ 


dC 


( ■'■Bq Vj) 2 _ 

= i{a 0 r Ul + B 0 Vj} io Q Uj + — U + 2i£{a 0 nij + 8 0 VjJWj 


- ^0 C5 

2{a 1 ru Q + v Q } -j { a o U 1 + 


3 0 V 1 


(4.14) 


whilst the z momentum and continuity equations give 


d? 


= 0 , 


,, „ , 6 0 V l) , d V lr *o c 
{a o u i + — rc'-x- 


- + e i v o 

a i u o + — — 


- 1C.J U_J 


le l V -1 


(4.15) 


so that Pq = C. It is important to point out at this stage that the terms 
proportional to U_^, V_^ in (4.14) are due to Coriolis effects; thus the 
structure of the neutral curve for stationary small wavenumber disturbances 
depends both on viscous and Coriolis effects. We can eliminate U^, 
from (4.14) and (4.15) to give 
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d 3 W, 



dW 

i{a 0 ru L + g 0 v^E , 2 + 2i5{a Q + B Q ^ }W ] _ 



- i ry o c ( - 

rUQ + v Q } j a 0 U 1 

B 0 ' 





~ U(0, /2s) 

0 U(0,0) ’ 


(4.16) 


where 



(4.16) 


We write the solution of (4.16) in the form 


W, = 


Yq _i/t | 9 

- i{a L ru Q + v Q } ^ - + A <Y 0 CF^s) 


8 


0 


2 i 60 

r —2-i 

1 + 1° 

Uq F 2 (s) 

r 

1 + -2 

. V 

U(0,0) j 


+ \ C , 


(4.17) 


where ki is a constant and Fi satisfies 


F'" - s F' + 2sF l = 1, F x (0) = F^oo) = 0, 


whilst F 2 satisfies a similar equation with the right-hand side of the 
differential equation replaced by U(0, /?s). In fact, it is straight- 
forward to express F^, F 2 in terms of integrals involving parabolic 
cylinder functions. It remains for us to satisfy = 0 at £ = 0; 

from (4.15) and (4.17) we can show that this condition leads to the 


eigenrelation: 
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2 10 O u 0 

Y 0 l 3 r 


1 + 


-2 

v 0 

=1 

u 0 


I, = 


iA ^ {a 1 ru„ + 


0 ' e i v oJ' 


Here the integrals I 3 , I 4 are given by 


X 3 = 


/ eu(o,e)de 

0 

2 u(o , 0 ) 


= . 599 , 


*4 = 


/ 0U (O,0)d0 

0 

U 2 (0,0) 


= .457, 


and (4.18) can be solved to give 


Y 0 = 


^ (‘ +^hJ l/2 r-‘ /2 - 

3 1 u o, 


-2 


a l _ 6 1 a 0 

8 o 


-2 

9 3/2/, ^ v 0 

2y o I 1 + =2 

u 0 

l? 0 -0 1 1/2 


r 1 / 4 


= 2.312 r 


-5/4 


We see at this stage that it is still not possible to find a 


1 


independently; however, it follows from above that <f> the angle 
the radius vector and the normal to the vortices is given by 


tan[ 7 r /2 -<)>] = 1.207 + 2.1312 e 2 r 1/4 + ,... , 


(4.18) 


(4.19) 


(4.20) 

and 0 ^ 
between 


(4.21) 
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1 '2 2 2 

whilst the total wavenumber — /a + 3 / r is given by 

e 


- l h 


1 / 2 , „2, 2 _ 1.224 r 

— /a + 3 /r - ^ + ,••• 

e e 


(4.22) 


The above expansion procedure can be continued in principle to any order 
and can take non-parallel effects into account in a self-consistent 
manner. We stress that (4.21), (4.22) have been obtained by taking the 

Coriolis effect into account; an Orr-Somraerf eld approximation to the full 
equations gives incorrect values for the second term in (4.21) and the 
first term m (4.22). The sensitivity of the structure of the lower branch 
modes to a combination of viscous and Coriolis forces means that, unlike 
the upper branch modes, for a more general three-dimensional boundary layer 
this class of modes might not even exist. Finally, we note that time- 
dependent modes with a sufficiently slow time scale are also possible and 
introduce a frequency into the eigenrelation (4.18). 


5* Conclusion 

The Reynolds number based on the boundary layer thickness, and the 

local azimuthal velocity of the disc is given by 


R 


A 



The inviscid modes have local wavenumber k. defined by 

A 


, / 2 2 , 2 

k^ = / a +g/r , 
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where the appropriate length scale is the boundary layer thickness. On the 
lower branch the local wavenumber k^ is defined by 


k 


A 


r V4 


/a 2 + 3 2 /r 2 


9 


so that (3.19), (4.22) are equivalent to 

k A = 1.16 - 14.4 R~ 1/3 + 


(5.1) 


and 


1.22 



+ 


) 


• • • 


> 


(5.2) 


respectively. Similarly (3.20), (4.21) become 

tan[ir/2 -<(.] = 4.26 + 29 R~ 1/3 + (5.3) 

and 

tan[ir/2 -<)>] = 1.21 + 2.31 R “ 1/4 + ,... . (5.4) 

Thus, if the neutral values are expressed in terms of R^, k^, and <p have 
no explicit dependence on the radial variable r. 

In Figures 2 and 3 we have compared the above asymptotic predictions 
with the numerical results of Malik (1985). The latter author solved the 
parallel flow approximation to (2.6) obtained by setting 3/9r = 10, 

3/30 = 16. Such an approximation is valid only for R -»■ ® but to the 
order shown in (5.1) - (5.4); our asymptotic results apply to the system 
solved by Malik. 
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We see in Figures 2 and 3 that there is satisfactory agreement between 
the asymptotic theory and Malik's results. Thus, the asymptotic approach 
will be a useful tool in finding the structure of the possible stationary 
modes in other three-dimensional boundary layers rather than having to 
solve the full parallel flow equations numerically. Similarly, the 
asymptotic theory could be used to identify the stationary modes which are 
likely to be important in a Navier-Stokes investigation of this problem. 

It is interesting to question why the lower branch modes have not been 
investigated earlier. The reason appears to be that in most experimental 
investigations of the disc problem, only the modes with <J> ~ 13° were 
observed. However, there is some discussion of modes with <f> ~ 20° in the 
paper by Federov, Plavnik, and Prokhorov (1976). These modes were found to 
exist closer to the center of the disc than the GSW modes and have a 
different vertical structure. Thus, it would appear that the lower branch 
modes perhaps bifurcate subcritically and therefore do not persist into the 
region where the GSW modes occur. Obviously, only a weakly nonlinear 
theory at least could settle this matter; however, it is interesting to 
note that Allen and Stuart (1985) have pointed out the possible existence 
of a subcntical mode with azimuthal wavenumber n = 2. 

The upper branch asymptotic results are again consistent with the 
results of Malik. The positive sign associated with the first correction 
term in (5.3) has important consequences. Thus, if (3.3) and (5.4) are to 
be connected at some finite Reynolds number the higher order terms in (5.3) 
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must be negative. This means that for some range of R the value of a 

A 

along the upper branch modes will stay close to the infinite value of about 
13°. This is exactly what Malik found numerically and even at the critical 
Reynolds number; <£ is still close to 13°. The wavenumber, however, changes 
much more along the upper branch; this presumably explains why GSW 
predicted $ so well but not the number of waves. 

Finally, we turn to the relevance of the lower branch modes m other 
boundary layer flows. At first sight we might think that our analysis is 
directly applicable to a two-dimensional boundary having zero shear stress 
at some position along the boundary. However, it is easily shown that the 
structure given in Section 3 is only applicable to boundary layers having a 
non-zero third normal derivative of the streamwise velocity component. 
This constraint effectively means that there are no neutral modes of the 
type found in Section 4 for spatially varying two-dimensional boundary 
layers. The Stokes layer velocity profile is another matter; at high 
Reynolds numbers the flow varies slowly in time, and the modes discussed m 
Section 4 are relevant to the times in a cycle when the shear stress 
instantaneously vanishes at the oscillating wall during the fluid motion. 

For three-dimensional boundary layers we expect that the lower branch 
modes are directly relevant. Moreover, it is of course possible that in 
such flows nonlinear effects might cause them to be more important than the 
GSW modes in the development of crossflow vortices. This matter can, of 
course, only be resolved by further calculations. 
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Figure 2. Comparison between the results of Malik (1985) and the 
asymptotic wavenumber predictions. 
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